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Abstract 

We consider the direct and indirect detection signatures of WIMPs x° i n kinematic 
regimes with a heavier, but nearly degenerate, charged state x^- F° r small splittings of 
O(10) MeV, the scattering of WIMPs off nuclei may be dominated by inelastic recombination 
processes mediated by the formation of (x~N) bound states, leading to a distinct signature 
for direct detection. These cross-sections are bound primarily by limits on the abundance of 
heavy isotopes, and may be considerably larger than the elastic scattering cross section in 
more conventional models. If the mass splitting is too large for recombination to occur, there 
may still be a significant resonant enhancement of loop-induced electromagnetic form-factors 
of the WIMP, which can enhance the elastic scattering cross-section. We also discuss how 
this regime affects the annihilation cross-section and indirect detection signatures, and note 
the possibility of a significant mono-energetic 7-signal, mediated by resonant processes near 
the ix + X~) bound state threshold. 
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1 Introduction 



The evidence for the existence of non-baryonic dark matter now comes from many sources 
and ranges over many distance scales pQ , from the rotation curves of galaxies, the dynamics 
of clusters, lensing data, and the characteristics of large-scale structure, to the features 
of the cosmic microwave background (CMB) fluctuation spectrum and the success of big 
bang nucleosynthesis (BBN). All of these pieces of astronomical data point to a similar 
cosmological density of dark matter, several times that of visible baryonic matter. Moreover, 
recent observations of the bullet cluster [2] have for the first time pointed more directly to 
the existence of dark matter through spatial separation of the baryonic and non-baryonic 
(or dark) components in the collision. 

This situation has been, and remains, one of the primary puzzles within particle physics 
and also one of the strongest motivations for physics beyond the Standard Model. It is 
remarkable that the theoretical shortcomings of the Standard Model with regard to the UV 
sensitivity of the Higgs mass point to new physics at or near the electroweak scale, and that 
a stable weakly interacting particle of this mass has just the right properties to be a dark 
matter candidate arising as a thermal relic from the Big Bang. There are also interesting links 
between various symmetries required to avoid large violations of approximate SM invariances, 
such as baryon number, and the stability required for such dark matter candidates. This 
apparent "naturalness" of weak-scale cold dark matter, or WIMPs, has led to a vast literature 
on the subject in models such as the super symmetric version of the Standard Model (MSSM) 
[lj, models with large extra dimensions [3], etc, and also the development of ground-based 
direct-detection facilities searching for the recoil of large nuclei from elastic scattering with 
WIMPs in the galactic halo U [5]. 

More recently, high precision data from WMAP |6J, and the negative results of the 
LEP and Tevatron searches for new physics at or near the electroweak scale has put some 
pressure on this seemingly natural scenario. In particular, in the MSSM with unification of 
the soft-breaking parameters at high energy scales (CMSSM), it is now well-known that a 
generic point in the parameter space, consistent with collider constraints, would typically 
lead to a relic density for the neutralino LSP which is orders of magnitude too high, while 
only a rather tuned region of parameter space remains in agreement with data [HE]. In 
fact, the simple freeze-out calculation of the neutralino WIMP relic density, which implies 
fl x h 2 ~ 3 x 1CT 37 cm 2 / (a ann v ) [Sj, leads to a value consistent with WMAP only in a region 
of the CMSSM parameter space which is already excluded by direct search bounds on the 
Higgs and chargino masses. Viable regions actually rely on additional, and to a certain 
extent accidental, features of the spectrum in order to enhance the neutralino annihilation 
cross-section. One such enhancement mechanism, the coannihilation [9] of the neutralino 
with a nearly degenerate charged slepton, allows the extension of the region of viable WIMP 
masses into the TeV energy domain. This issue of tuning is of course not unique to the 
relic density of dark matter and is now a more generic problem for models of new weak 
scale physics following the absence of new discoveries at LEPII for example. Indeed, the 
relic density requirement and the apparent absence of new physics right at the weak scale 
suggests a prominent role for mechanisms which allow for enhanced annihilation of heavier 
WIMPs in coming years. 
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With this motivation in mind, in this paper we would like to focus on one particular 
enhancement mechanism, namely coannihilation with a nearly degenerate charged state, 
and explore how extreme limits of this kinematic regime may alter the predictions for direct 
and indirect searches. We will refer to this kinematic scenario as a pseudo- degenerate WIMP, 
having in mind that the WIMP Xi will be separated by a small mass gap of Am ~ 0(1 — 100) 
MeV from an excited state X2 or X2 ■ For a WIMP mass of (9(100 GeV) this splitting is 
actually much smaller than is required for coannihilation, for which a relative splitting of less 
than 5% is sufficient. Whether such near-degeneracies are natural or not clearly depends on 
the precise nature of the model at hand, but it is a particularly interesting kinematic regime 
in the case that \i is charged for the following reason. For Am < (9(20 MeV), inelastic 
scattering of Xi with heavy nuclei becomes possible in which capture of X2 occurs to form 
the bound state [Nxz)-, and this process can dominate the cross-section if xi behaves in other 
respects as a conventional WIMP with a suppressed elastic scattering cross-section. This 
process, since it involves charged exchanges, can lead to significantly different signatures 
in direct detection experiments. If Am is larger than about 20 MeV, this recombination 
process is no longer possible, but there can still be a significant resonant enhancement of 
the elastic scattering cross-section. Such enhancements of elastic scattering are in some 
respects the trade-off for going to a kinematic regime where the annihilation cross-section 
is similarly enhanced. For very small splittings, the annihilation cross-section itself may 
also be dominated by resonant processes, since it is near threshold for the formation of the 
(x + X~) bound state, which in turn can result in an enhanced production of 7-rays. This is 
particularly so in the galactic environment where the characteristic velocities are quite low. 

In what follows we will explore these issues in turn. In Section 2, we introduce the 
general kinematic scenario and some classes of pseudo-degenerate models. In Section 3, 
we consider first the recombination with nuclei and compute the simplest charged current 
and electromagnetic captures to the ground state, and comment on various constraints from 
e.g. terrestrial heavy isotope searches. We then turn to the generic resonant enhancement 
of the elastic cross-section. In Section 4, we consider the impact of resonant processes in 
annihilation, and we conclude in Section 5 with some additional remarks on alternative 
detection signatures. 

2 Pseudo-degenerate WIMPs 

Following the motivation outlined in the preceding section, we will consider generic WIMP 
scenarios in which the dark sector has some substructure, in the form of at least one excited 
state X2 nearly-degenerate with the WIMP xi which we also take to have a mass of O(0.1 — 
1 TeV), with the lower limit imposed by searches at LEPII. A near-degenerate neutral excited 
state was considered previously for different reasons in [10], but for efficient coannihilation 
X2 should have stronger interactions with the SM than xi an d for the present paper we 
will generically assume that X2 is electromagnetically charged, X21 although a neutral state 
charged under SU(3) C would also fall into a similar class. These simple requirements define 
the class of pseudo-degenerate WIMPs to be studied in this paper. 
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2.1 Models 



To provide a somewhat finer classification of pseudo-degenerate scenarios, we will generally 
consider two classes of models, determined by the dominant interactions with gauge bosons. 
The WIMP xi, whether fermion or scalar, is required to possess no diagonal tree-level spin- 
independent coupling to SM gauge bosons due to constraints on direct detection. Even a 
very small, O(10~ 2 ), spin-independent coupling to the Z is still sufficient to produce an 
elastic scattering cross-section with nuclei in excess of the current bound. The nullification 
of couplings to Z can be achieved either by a careful charge assignment in the WIMP sector, 
or by requiring that xi an d its charge conjugate field x\ are the same. In other words, in 
order to suppress the spin-independent scattering cross section, xi should be a real scalar if 
WIMPs have spin 0, or a Majorana fermion if xi has spin 1/2. We have defined our scenarios 
such that xi is electromagnetically charged, and thus the primary distinction we can place 
on the xi ~~ Xi WIMP sector is on the type of charged current between these two states. 
Restricting our discussion to tree-level couplings, we introduce two model classes. 

• Type A — The first scenario we will consider will allow for a charged current interaction 
with W-bosons. In other words, the pseudo-degenerate WIMP sector possesses an off- 
diagonal vector current J£ lX2 so that, 

Ant = J&,W+ + h.C, (1) 

and consequently Xi an d Xi are either both bosons or fermions. A familiar MSSM 
model in this class would involve a neutralino xi with a near-degenerate chargino xf- 

• Type B — The second scenario will be defined by the condition that xi an d Xi are 
respectively a scalar and a fermion or vice versa. In this case the current between xi 
and xi is fermionic and couples to the SM charged leptons ip 

c int = r x ; X2 ^ + h. c .. (2) 

A characteristic MSSM example in this case would be a neutralino Xi with a near- 
degenerate stau Xi- 

In the remainder of this section, we outline the relevant kinematic regimes of interest. 
2.2 WIMP-nucleus binding energies 

Coannihilation does not impose undue levels of tuning on the mass spectrum; a splitting 
between the WIMP and excited states of 0(5%) is generally sufficient. This does not signifi- 
cantly alter the expectations for interactions with baryonic matter, and we will be interested 
in a more extreme limit in which the splitting is in the MeV range, and thus comparable to 
the Coulomb binding energy of X2 with nuclei. We will defer any discussion of theoretical 
motivations for this pseudo-degeneracy, but as a benchmark point to normalize the kinematic 
regimes to be discussed below, it is worth noting that if Xi an d xf were truly degenerate 
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54 


21.7 


18.4 



Table 1: Estimates for the binding energies of the state (A r x 2 ) assuming a gaussian and step-like nuclear 
charge distribution for several relevant elements. 



in the dark sector through some symmetry, then the interaction of xf with the SM gauge 
bosons would naturally imply a splitting of 0(100 MeV) [HjQ 

Denoting the mass difference as Am = m X2 — m Xl , we explore the upper bound on Am 
that leads to the formation of stable (Nxz ) bound states. In the limit m x 3> m^, the 
binding energy depends only on the nuclear mass, and is naively given by the analogue of 
the Rydberg energy, —E b ~ {Za) 2 m^ /2. More precisely, since the characteristic radius of 
the orbit ~ [Zamfq)~ x generally lies well within the nucleus, we obtain a better estimate 
by solving the Schrodinger equation with a given charge distribution inside the nucleus. This 
leads to the results shown in Table 1 for several elements that will be relevant in this paper. 
Two types of charge distribution with the same (rf), Gaussian and step-like, are employed 
that can be viewed as two extreme approximations of a more realistic nuclear charge density. 

Table 1 reveals three distinct kinematic regimes: 

[(I) : recombination with light elements] < Am < 4 MeV, (3) 

[(II) : recombination with heavy elements] 4 MeV < Am < 20 MeV, (4) 
[(III) : resonant loop enhancement] 20 MeV < Am < 100 MeV. (5) 

In regime (I), stable bound states with light nuclei up to oxygen may be formed. Some 
models in this kinematic regime are significantly constrained by searches for anomalously 
heavy isotopes of carbon and other light elements as discussed below. Regime (II) is equally 
interesting, as it opens the possibility of 'recombination' processes where WIMPs may com- 
bine with the heavy elements inside the detectors used for the direct searches of dark matter. 
Finally, kinematic regime (III) arises when Am is too large for recombination to occur in 
inelastic scattering, but still small relative to the WIMP mass, and the elastic scattering 
cross-section can be resonantly enhanced through loop processes, with an off-shell Xo. a P" 

1 It is also worth noting that if \i an d X2 arise from the low-energy hadronization of a single relic state 
charged under SU(3) C , e.g. as hadronized squarks [12], they can naturally be split by O(MeV) as members 
of an isospin doublet. In such cases, recombination with nuclei could occur through strong rather than 
electromagnetic interactions, but this may be subject to more stringent constraints than the scenarios for 
electromagnetic capture that we consider here. 
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Figure 1 : The weak charged-current capture of Xi to form the bound state (N\2 ) ■ 

pearing as an intermediate (virtual) state in scattering. The upper limit in this case is not 
strict, but simply gives an estimate for the level beyond which these processes are less rele- 
vant; 100 MeV is the characteristic centre-of-mass momentum for interactions of a 100 GeV 
halo WIMP with terrestrial nuclei. 

3 Pseudo-degenerate WIMP- nucleus scattering 
3.1 WIMP-nucleus recombination 

One of the more interesting processes that becomes kinematically accessible for the mass 
splitting in (j3j), is the 'recombination' of WIMPs with nuclei, and indeed such inelastic 
scattering can dominate the cross-section if xi is otherwise relatively inert, as would be 
characteristic of a WIMP. In this section, we explore several classes of interactions. 

3.1.1 Weak charged-current recombination 

For models of type A, a natural class of capture processes will proceed via weak currents, 
the simplest example of which is shown in Fig. [TJ Xi + n ~ * X2 + P- Depending on the 
value of Am, the recombination process may occur directly to the ground state of the bound 
WIMP-nucleus system, or to an excited intermediate state that will subsequently decay 
to the ground state through emission of 7's and/or neutrons depending on the nucleus in 
question, 

x o + N (z) _ (tf(*+i) x -)« _ {N^ X ~) + ( 7 , n, . . .). (6) 

This process may receive a resonant enhancement but of course is also subject to the details 
of nuclear binding and so may not be energetically allowed for many light nuclei, depending 
on the relative binding energies of N^ z ' and N^ z+1 \ For example, the capture ([6]) is not 
possible for such abundant nuclei as 12 C, 14 N, and 16 0, because the (Z + 1, A) nuclei are too 
massive. In addition to (jB]), there always exists a /3 + -type process, 

x o + N (z)^ {N (z) x - ) + e+ + Ui (7) 

which is clearly non-resonant. 

The non-resonant contribution ((7|) is relatively easy to estimate, and to refine to a full 
calculation if needed. For our purposes it suffices to estimate the non-resonant capture cross 
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section by considering the overlap between the scattering and bound state wavefunctions, 
i.e. 

^non-res * | (^«* |^> ? X Fq, (8) 

where T Xl ^ X2 is the "decay" width of xi to xi with the energy release Q e s = \Eb\ — m e — 
Am. Neglecting the Gamow factor Fq for now, saturating the wave-function overlap by 
characteristic nuclear scales, and taking T X1 ^ X2 ~ 10 _3 Hz(<5 e fr/lMeV) 5 , we immediately 
discover that the rate for this process is extremely slow, 

Wn.« ~ l(T 3 Hz x (2 fm) 3 x ( * ~ lO^cm 2 ( * , (9) 



^lMeVy V 1MeV , 

even for the largest values of Q c s- 

Computing the cross-section for process (jHJ) with an intermediate excited nuclear state is 
rather nontrivial. However, to get a reasonable estimate, it is useful to think of the excited 
state {NX2)* as a resonance, for which Fig. [1] characterizes the entrance width rV- Since Tw 
relates to a weak process, it is much smaller than the decay width induced by electromagnetic 
and strong interactions, and we can use the Breit-Wigner formula to write, 



R 



where g* refers to the angular momentum multiplicity factor, q cm is the momentum in the 
WIMP-nucleus centre-of-mass frame, and Er denotes the resonant levels of (iVxj) which 
may be formed in the capture process. 

The spread in kinetic energy of the colliding WIMP-nucleus system is on the order of 
M R (N, x) v2 i where v is the relative velocity and Mr the reduced mass. Thus, for nuclear 
masses in the interval 10 - 100 GeV, the typical spread in kinetic energy is 10 - 100 keV, 
which is somewhat smaller than the characteristic spacing of bound state energy levels of 1 
MeV. This level spacing becomes even denser due to the additional Coulomb excitations in 
the (NX2) system, and for the purpose of obtaining an estimate it is reasonable to assume 
that AE ~ 0.5 MeV. The decay width will vary quite significantly depending on Am and the 
nuclei participating in capture process. For electromagnetic decay widths in bound states 
with light nuclei it is reasonable to expect T ~ 0(1 eV), while a neutron decay width in the 
capture of a WIMP with small Am by a large nucleus can easily reach O(100 keV), and may 
thus be quite comparable to the spacing of energy levels. For the latter case, we approximate 
the capture rate by assuming a characteristic energy denominator scale in (fTOl) on the order 
of AE, 

VvesV ~ -TJ2— X ■ (11) 

M R v 500 keV 

To complete the estimate we need to evaluate the weak entrance width. Its dependence 
on the main parameters in the problem is captured by the following scaling, 

GlMlv 
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Figure 2: Another charged current capture process of Xi forming the bound state (-/Vx 2 ) and radiating a 
neutrino. 



where tq is the characteristic distance scale entering in the evaluation of the matrix element, 
that for the purposes of this estimate we take to be of order the nuclear radius of 4 fm, 
and we have also taken the coupling of X1X2 to W to be 0(1). The estimate (fT2|) may be 
viewed as the maximum possible weak rate, as it supposes a high degree of coherence among 
nucleons. Indeed, the M\ dependence suggests a coherent response of all the nucleons at 
characteristic momentum transfers of order Mrv, which might lead to an overestimate for 
large nuclei. One finds that a parton-level calculation of the width, which effectively ignores 
coherence effects completely, leads to a result that is about an order of magnitude smaller. 

Combining the two pieces together, (TTTj) and ffl2l . we arrive at an estimate of the capture 
rate given by 

a„v ~ 3 ,^, v , ~ 1CT 38 cm 2 . (13) 
Tq x (500 keV) 

For light nuclei where the outgoing electromagnetic width is on the order of 1 eV, the 
rate is expected to suffer from a further suppression by T em /AE, i.e. another five to six 
orders of magnitude. We observe that the result is similar to the natural scale for this 
weak process, and thus does not appear to be significantly enhanced over more conventional 
elastic scattering of Xi- However, we emphasize that this treatment of the nuclear aspects of 
the calculation has been very cursory, and a more complete treatment may lead to further 
enhancement (or suppression) factors. 

The capture process above is quite distinct from a characteristic nuclear recoil event due 
to elastic WIMP scattering because of the large O(MeV) scale energy release during the rapid 
decay of the excited bound state. Naively this makes it very difficult to observe in existing 
direct detection experiments which have detectors tuned to detect a fiducial recoil energy of 
no more than 100 keV. For this reason, it is interesting to consider inelastic processes that 
would have a signature somewhat closer to a standard recoil event, and one possibility is 
shown in Fig. [5J for the special case in which \2 is a scalar partner of the r as in the MSSM. 
The reverse process, a nuclear-assisted stau to neutralino decay was discussed previously in 
[T3] . Such a capture process is analogous to ([6]), with the important distinction that the 
bound state de-excites through radiating neutrinos which would escape the detector leaving 
just the recoil signal to be observed. However, this process is further suppressed and may 
not significantly impact the cross-section purely through elastic scattering of Xi- 
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Figure 3: The electromagnetic capture of Xi to form the bound state {N\ 2 ) an d radiating a positron. 
3.1.2 Electromagnetic recombination 

For models of type B, which lack a weak charged current, alternative recombination channels 
are open if the charged states represent bosons and the neutral states fermions, or vice versa. 
In particular, the electromagnetic processes shown in Fig. |3] become possible. In the early 
universe the abundance of charged states is then rapidly depleted via processes such as 
xf ~ ► X? + e± i while in the current epoch \~2 ma y De regenerated via recombination in the 
form, 

Xi + N^(Nxz) + e+, (14) 

which may occur given the appropriate kinematics (l3f4l) . where N is again a generic nucleus, 
and (N\2 ) a stable bound state. Process (|T4"1) may appear to violate lepton number, and 
thus be suppressed. However, this need not be the case if either %t or Xi carry lepton 
number as happens e.g. in models where Xi is a sneutrino, and Xi a chargino, or when xi is 
a neutralino, and X2 a charged slepton. In the widely discussed neutralino-stau scenario, the 
coupling to the positron would actually represent a flavour- changing process in the lepton 
sector, which may be somewhat suppressed but does not have to be vanishingly small. 

Process (fT4l) is an interesting variant of standard radiative recombination of a nucleus 
with an electron, in that Coulomb interactions are present in the final rather than the 
initial state, and since the positron has to tunnel out of the nucleus the cross-section will be 
Gamow-suppressed, particularly when the positron is non-relativistic. To be concrete, we 
will consider the situation where Xi is a scalar and Xi a Dirac fermion, with the effective 
vertex, 

£ = 9XieX2 + h-c. (15) 

The free decay width T X2 of X2 in the early universe is then sufficiently fast to avoid problems 
with BBN provided g 2 > 10" 17 [H]. 

We estimate the capture rate by considering the overlap between the scattering and 
bound state wavefunctions as in ([8]), i.e. a ~ r xi _ >X2 |(^ at |?/;^)| 2 x F G where T X1 ^ X2 is the 
free decay width (for Am < 0), and Fq is the Gamow factor associated with the interaction 
of the radiated positron with the nuclear potential. In the centre-of-mass frame of Xi and 
N, the positron energy is E e + m \Ef,\ + l-E^ml — Am + m e ~ 5-15 MeV, while the nuclear 
barrier is characteristically of 0(Za/R^) that may reach 20 MeV for large nuclei. In that 



case the positron will have to tunnel out of the Coulomb barrier leading to a suppression in 
the cross-section, i.e. F G <C 1. We will discuss this point further below. 

In what follows we consider the capture to the ground state that maximizes the kinetic 
energy of the outgoing positron and thus minimizes the Gamow factors. The bound state 
wavefunction can be taken as that of a 3D harmonic oscillator, with the inner non-Coulombic 
part of the nuclear potential given by V(r) = (Za) / (2Rn) ((r / Rn) 2 — 3). Computing the 
overlap and converting to nuclear parameters, we obtain 

g 2 {E e + + m e )\p e +\ 3 / a B \ 3/4 3/2 / n 2d 3/2 iM xx „ , 1fi s 
{(TrecV) = — R N J 8^ exp {-(jiv) J x F G , (16) 

where fi = M B (N,x) is the reduced mass of N and xi, = l/(Zafi) is the 'Bohr radius', 
and Rn the nuclear radius. 

If we suppress the detailed dependence on the bound state wavefunction, then in the 
limit of a relativistic positron and assuming a B ~ R N , this reduces to 

{W) ~ gfa%^, (17) 

which exhibits the inverse scaling of the inelastic cross-section with the WIMP momentum. 
The most significant numerical correction to ffTTj) arises from the fact that in practice a B <C 
Rn for the nuclei of interest. The exponential factor in fflBI) . associated with the ground 
state wavefunction, provides a correction of 0(20%) for large nuclei. Evaluating the estimate 
(JTTj) for WIMP capture on light nuclei, such as carbon or nitrogen, we take a B ~ 1 fm, and 



E e + ~ 1 MeV, to obtain the recombination rate 



i-\ x (il eV 



WvccV) light nuclei ~ 10 34 Cm 2 X ( ) X ( - ~ ' ) . ( 1 > 



For heavier nuclei, the Gamow factor F G can no longer be ignored and leads to significant 
suppression if the positron energy is well below the nuclear potential barrier. The present 
situation differs slightly from conventional non-relativistic tunneling, where F G ~ eT G with 
G = j c dr\k\ ~ 2nZa/v, since for generic values of Am the positron will be relativistic even 
for energies well below the barrier. Consequently, Za/v will not be too large and we can 
ignore Coulomb resummation effects after tunneling has occurred and simply account for the 
tunneling factor. This we estimate as G rel = Re j c dr\kb\ using the relativistic momentum k. 



>i, 



under the (wide) nuclear potential, which leads to G rel ~ 2irZaE/k provided the positron 
is only mildly relativistic. The latter constraint is necessary to be in the "wide barrier" 
regime which restricts the positron to have 7 < 0(5). For the values of Z relevant here, 
the cross-section is still somewhat suppressed for generic 7-factors of the emitted positron, 
but far less so than in the non-relativistic regime. The impact is most significant for our 
maximal choice of Am as the positron is then barely relativistic even with capture to the 
ground state, and Coulomb effects may also become important. For a radiated positron with 
7 ~ 1 - 2, we find F G ~ 1(T 3 - 10 _1 . 

Inserting characteristic numerical scales for heavy nuclei into (TIBl) . taking E e + ~ 10 MeV 
and setting F G ~ 1CT 2 , we can convert this into an effective cross-section per nucleon, the 
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(19) 



where Mr denotes the reduced mass, and which for g 2 /(47r) ~ a w is interestingly a few 
orders of magnitude above the current direct detection constraint on elastic scattering, e.g. 
with Ge or Xe HE]. 

Its also worth noting that an additional inelastic channel is triple recombination with the 
use of a i^-shell electron, 



This process does not require tunneling through the potential barrier, but is suppressed by 
the small probability of finding the i^-shell electron near the nucleus, {Zam e R^Y . This 
channel is expected to remain subdominant except for the case of very small E e +, i.e. with 
Am very close to the capture threshold. 

3.1.3 Constraints from BBN and heavy isotope searches 

As is well-known in the context of charged dark matter, there are stringent constraints on the 
presence of long-lived charged particles surviving from the Big Bang. It has recently been 
recognized that charged states with lifetimes longer than (9(10 3 s) can also severely disrupt 
the successful predictions of BBN through the catalysis of interactions that, for example, 
would overproduce 6 Li and 9 Be by many orders of magnitude [TJJ [15]. For this reason, we 
will simply assume here that the state xf has a sufficiently small lifetime to decay before 
nucleosynthesis: 



This can be a significant constraint on realizations of this scenario within the MSSM, where 
for example X2 = in the absence of a flavour-changing stau-neutralino-electron coupling, 
the condition (|2"T|) is equivalent to requiring that the mass splitting between stau and neu- 
tralino is in excess of 100 MeV. However, even an extremely small coupling to the electron, 
g ~ O(10 -8 ) would allow (T2~Tj) to be satisfied without imposing significant constraints on 
Am. 

More direct constraints on the mass splitting Am then arise from the terrestrial limits on 
the abundance of anomalous heavy isotopes. Here the strongest constraints apply to light 
elements, which are in some sense also the most reliable, arising from well-mixed gaseous or 
liquid media. The most stringent is the abundance constraint on heavy hydrogen (He^j), 
in the form of heavy water [16], which however is only relevant here if the mass splitting is 
very small, less than 350 keV. There are also quite stringent constraints on (C%j) (relative 
abundance /# < 10~ 14 per nucleon) and (N%j) (fc < 10~ 20 per nucleon) [T7] . 

To make use of these constraints, we connect the abundance of the parent nuclei that 
undergoes the capture process to that of the anomalously heavy daughter nucleus, a bound 
state with X2, 



X 1 + e- + N-+ (NxzT - (JVxa ) + 7- 



(20) 



r X2 <0(l)s. 



(21) 




(cr rec u)n DM T f 



(22) 
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Here r exp is the exposure time, and 7Zdm is some average dark matter density near the solar 
system, which we shall assume to be equal to the local WIMP number density, tzdm ~ 
3 x (TeV/m x ) x 10~ 4 cm -3 . Specializing to the case of electromagnetic capture (1141) on 
nitrogen, and using the limit on the anomalous carbon abundance, fc < 10~ 20 per nucleon, 
we arrive at the following parameter constraint, 



In deriving this constraint we also assumed that carbon and nitrogen are almost equally 
abundant in nature, which is rather conservative: the carbon used in Ref. [17] presum- 
ably spent the majority of its time in the atmosphere, where nitrogen is significantly more 
abundant. Nonetheless, the constraint (1231) is very significant, and for example implies the 
absence of a large coupling to the first generation of leptons in the neutralino-stau model. 
It clearly disfavors models of type B with mass splittings below ~ 3 MeV, when the binding 
to nitrogen becomes energetically possible. 

It is remarkable that the abundance constraints on anomalous isotopes of heavier el- 
ements are considerably weaker, while there is also greater uncertainty over whether the 
samples tested have a characteristic exposure time. Even interpreted in the most conser- 
vative manner, existing bounds e.g. on heavy isotopes of Au [18], from (Hgx^), lead to 
relatively mild constraints that are well below the estimates for the cross-sections obtained 
in this section. For this reason, the range quoted in (j3J) seems perfectly viable, although it 
would of course be interesting were further searches to be performed for heavy isotopes of 
elements with Z > 10. 

Thus, we conclude that for models of Type B heavy isotope searches do impose a sig- 
nificant constraint but only in connection to binding with light elements and do not overly 
constrain the cross sections for Am in the range (jl]). For models of Type A, where only 
weak charged currents are allowed, we noted earlier that the nuclear binding energies do not 
permit capture by the light C, N, and O nuclei. Thus, since the capture rates are too slow 
to provide a significant constraint from searches for anomalous isotopes of heavy nuclei, the 
entire range of Am including the smallest values down to an MeV is allowed. We note in 
passing that models with weak currents, Am ~ few MeV, and the X2 lifetimes of ~ 2000 
seconds alleviate the known problems with the primordial abundance of lithium |15j . 

3.1.4 Constraints from direct detection and from annihilation in the sun 

The existing limits on the elastic cross-section per nucleon have reached impressive levels of 
order 10 -43 cm 2 for a 100 GeV WIMP and further progress is anticipated in coming years. 
However, its worth bearing in mind that the detection strategy in this case relies crucially 
on a characteristic recoil signal with energy of order 50 keV. In this context, its interesting 
to consider if direct terrestrial dark matter searches are indeed sensitive to electromagnetic 
recombination processes of the type shown in (fl4l) . At first sight, it is unclear if the recoil 
during the capture process would be observed, due to the additional MeV-scale energy release 
associated with the subsequent positron annihilation in the detector. For this reason, we 
cannot immediately place constraints on the size of the coupling g from this source. Indeed, 




for Type B models with Am < 3 MeV. 



(23) 
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Figure 4: A schematic representation of the resonant effects we will consider, where the shaded region 
denotes the interaction with the nucleus. 

for the majority of detectors, events with such a large energy release would be identified with 
background and rejected. 

Besides the direct detection strategy that is currently focused on the search for WIMPs 
elastically scattering off nuclei, another traditional method to look for WIMPs is by searching 
for highly energetic neutrinos coming from WIMP annihilation in the solar core. The density 
of WIMPs in the solar core can be many orders of magnitude larger than rirjM due to 
gravitational capture through down-scattering of WIMPs on nuclei inside the sun, with 
an efficiency again controlled by the elastic scattering cross section. However, the pseudo- 
degenerate WIMP regime can significantly modify the sensitivity of this method. Indeed, if 
Am is less than about 10 MeV, sufficiently small to allow recombination with iron, a trapped 
WIMP will preferentially undergo iron capture rather than participate in direct annihilation 
with another WIMP. This stems from the fact that the density of iron in the solar core is 
many orders of magnitude larger than the density of WIMPs, and that the capture rate, (1121) 
or ffTTj) . is quite comparable to or even larger than the annihilation rate. Simple estimates 
show that for a typical capture rate of (av) ~ 10 _40 cm 2 , a trapped WIMP recombines with 
an iron nucleus within 10 8 seconds which is a much shorter time scale than for annihilation 
with another WIMP. Therefore, most of the WIMPs inside the sun would be in bound states 
with iron and therefore shielded from annihilation by a large Coulomb barrier, which would 
essentially cut off the flux of ultra-energetic neutrinos. 

3.2 Resonant enhancement of EM form-factors 

The characteristic momentum for an (9(100 GeV) WIMP in the galactic halo scattering with 
nuclei is m x v ~ 100 MeV, which implies that even if the mass splitting Am between xi an d 
X2 is too large to allow for recombination with nuclei, significant resonant enhancements 
from the off-shell \2 are possible in the elastic scattering cross- sect ion. This is the scenario 
we will explore in this section, having in mind a splitting in the range so that \2 appears 
off-shell with virtuality q 2 ~ mAm in the scattering process as shown schematically in Fig. 4. 

The precise nature of the amplitude will depend on the mechanism via which the WIMP 
interacts with the nucleons. The general XiN ~ > XiN elastic scattering amplitude will 
involve many contributions, but since the WIMPs in the galactic halo are non-relativistic, 
one can summarize the dominant contributions through the leading low-g 2 components of 
form-factors of the WIMP with various currents to which the nucleon couples. Imposing 



12 



minimal assumptions on xi, namely the absence of tree-level vector couplings to 7, Z, and 
gluons, there are a number of form-factors which may play a role in interactions with the 
nuclei in the detector. To keep the discussion manageable, we will restrict our attention to 
electromagnetic form-factors, although one should bear in mind that much of the discussion 
will translate directly to couplings to other SM gauge bosons, and in various scenarios gluonic 
or weak form-factors may be of more importance. It is also important to keep in mind that 
these form-factors are often only well-defined off-shell and thus in practice this is just a 
convenient means of keeping track of the leading parts of the scattering amplitude in an 
appropriate kinematic regime. 

Recall that for a neutral state x °f arbitrary (nonzero) spin, if we impose T-invariance, 
we can expand the matrix element of the electromagnetic current [19] as follows in powers 
of the momentum transfer q = p — p', 



-\Q(Sq)\ + 



Xiip), (24) 



where is the spinorial generator in the appropriate representation (i.e. for a spin-1/2 
particle) and the corresponding spin 4- vector is = ie tlup<T T,p a q 1 , / (4m x ) . The small-g 2 limit 
allows an expansion of the general form-factors and the terms on the first line contain the 
leading moments for a spin-1/2 particle - the charge radius rf>, the anapole moment a which 
is P-odd, and the magnetic dipole moment /1. For spin S > 1, there are also higher order 
multipoles and the second line contains the leading term in this limit, namely the electric 
quadrupole moment Q. Many of these moments vanish given specific constraints on xi, while 
on the other hand if we relax the constraint on T-invariance, there are in addition the electric 
dipole moment, the magnetic quadrupole moment, etc. but these will not be needed in what 
follows. However, particularly in the case of scalar particles where (with the exception of 
an appropriately defined charge radius) these moments do not exist, we need to go beyond 
one-photon exchange and consider two-photon processes. We will again consider the small 
q 2 limit, in which case the leading constant parts of the form-factors are the electric and 
magnetic susceptibilities T]e and t]b, and write the matrix element in position space, 



£ 27 = xi^Tf 



X + • • • (25) 



Once again, if we allow for T-violation, there are also mixed polarizabilities given by x^Fx, 
but we will drop these contributions in what follows. 

At the non-relativistic level, the relevant constant coefficients in the form-factors are 
conveniently assembled into a Hamiltonian describing the interaction of the non-relativistic 



2 For particles with no vector coupling to the photon, the anapole moment a (and also the charge radius r 2 D ) 
are strictly defined in a gauge invariant manner only through an external current. They determine a coupling 
to an off-shell photon, as is relevant to nucleon scattering, and one should also include the corresponding 
off-shell Z coupling for gauge invariance. We will ignore this issue as in practice for the momentum scales 
relevant here this effect is negligible. 
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Figure 5: A contribution to the scattering cross-section of \i with a nucleon, mediated by the multipolc 
moments (on the left) or susceptibilities (on the right). 

WIMP with a slowly varying electromagnetic field. The general form of this Hamiltonian 
was considered previously in [20], and with the restriction to T-invariance as above, takes 
the form, 

Hl even = -fJ,B-S-ai-S-~Q ij d i E j 

-\er 2 D V • E - l - VE E 2 - l -r, B B 2 + ■■■ (26) 

with Qij = Q(S(iSj) — (2/3)SijS(S + 1))/(S(2S — 1)), in terms of the various multipole 
moments and susceptibilities above. The moments on the first line apply only for nonzero 
spin. If the particle is a Majorana fermion then furthermore Qij = fi = r 2 D = 0. The 
vanishing of the effective charge radius of a WIMP has significant consequences for WIMP- 
nucleus scattering, as r 2 D contributes to spin-independent scattering enhanced by a very large 
nucleon coherence factor. 



3.2.1 Resonant enhancement of loop processes 

The cross-sections for non-relativistic scattering due to these moments are listed for example 
in [20]. The main point we wish to emphasize here is that even if these moments are zero at 
tree-level, for pseudo-degenerate WIMPs they naturally arise at loop-level and indeed may 
be resonantly enhanced. The characteristic diagrams are shown in Fig. 

Let us focus first on models of type B, where the WIMP xi is a fermion with a nearly 
degenerate scalar partner \2 which appears resonantly in the scattering process. Its con- 
venient to denote the "moments" generically as M^) where d is the mass dimension of the 
corresponding operator, i.e. d — 1 for /i, d = 2 for r 2 D and a, while d = 3 for t)e- An inspec- 
tion of the loop diagrams indicates that the decoupling with the WIMP mass m x is generic, 
i.e. Mfj\ ~ l/ m x' gi ven that the extra internal fermion state is light, consistent with reso- 
nance. This scaling can be straightforwardly understood by taking the non-relativistic limit 
for the heavy fermions. This leads to a factor of v /m x at each vertex, while the heavy scalar 
propagator leads to a factor of 1/tu 2 . Moreover, in the resonant limit the loop is infra-red 
divergent and thus the result is enhanced by the small mass scale Am which sets the cutoff 
momentum scale, provided a light SM fermion with m < Am is available. Otherwise, one 
must replace Am with the corresponding mass. E.g. in the case of a near neutralino-stau 
degeneracy, the infra-red cutoff is fixed by m T . Thus in the generic case, the moments scale 



14 



in the following way in the resonant limit, 




(loop factor) 



x 



1 



(27) 



m 



(Am) 



X 



where in fact the loop factor is also enhanced by tt due to the infrared divergence. This 
formula is quoted for zero momentum transfer, and in practice for q 2 > m x Am > (m x v) 2 the 
EM form-factors will decrease as a function of q 2 . This scaling suggests that the enhancement 
is all the more impressive for the susceptibilities, but one must bear in mind that these 
amplitudes are 0(a 2 ) and more importantly involve a chirality flip and so require a double 
degeneracy of xi with both \2L and \2R- 

We have studied the contributions of these enhanced amplitudes to the scattering cross- 
section and, although significant, the results in most cases do not reach the level of the 
current direct search bounds. For this reason, we will not present full formulae, but simply 
comment on a few interesting cases. 

• Charge radius — The exception to this general conclusion involves the charge radius 
in any scenarios in which it is present at 1-loop, arising e.g. from the diagram on 
the left of Fig. [5j A straightforward calculation leads to a result of the form r 2 D ~ 
a/(m x Am) where the order-one normalization depends on the precise field content. 
For Am ~ (9(100 MeV) this leads to a (per nucleon) elastic scattering cross-section 
[20] on 72 Ge for example of order 10~ 38 cm 2 which is several orders of magnitude above 
the current bound, and close to the tree-level Z-mediated cross-section. Therefore, 
this can provide a rather stringent constraint on scenarios which, for example, avoid a 
coupling to the Z by choice of the SU(2)l representation. 

• Anapole moment — For charge-less Majorana fermions, there is a single nonzero 
multipole moment - the anapole moment [21] - and in the scenarios considered here 
it can be quite sizable, consistent with the above result for the charge radius. As an 
explicit example, one may again have in mind a near-degenerate neutralino and stau in 
the MSSM. However, as the anapole is an axial moment, in elastic scattering it couples 
to the spatial nucleon current and so the cross-section is suppressed by 0{v 2 ) ~ 10~ 6 
[20J which pushes it somewhat below the current level of sensitivity. 

• (Chromo) electric polarizability — The velocity suppression which afflicts the cross- 
section derived from the anapole moment does not apply to the polarizabilities. How- 
ever, these operators require a chirality flip and so for resonant enhancement demand 
a double-degeneracy of the left- and right-handed \2 states with \i- F° r weakly inter- 
acting states, e.g. ft and tr approximately degenerate with a neutralino (bino), the 
diagram on the right of Fig. pleads to xe ~ 10~ 2 ct 2 /(m x (Am) 2 ) in accord with the 
scaling of f[2~Tj) . However, even in this degenerate regime which resonantly enhances the 
amplitude, the cross-section for nucleon scattering is negligible due to the small electric 
energy density in the nucleus (N\E 2 \N) ~ 20 MeV. However, there is an analogous 
chromoelectric polarizability, rjc, replacing the photons in Fig. [5] with gluons which is 
far more significant since (N\a s G 2 \N) ~ m^. The corresponding contribution to the 
scattering cross-section in the resonant regime can be as large as <J( Pt n) ijlc) ~ 10~ 38 cm 2 
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Figure 6: A resonant contribution to the annihilation cross-section of X1X1 - > 27, mediated by the (X2X2 ) 
bound state. 

for Am ~ (9(100 MeV). This result is contained as an appropriate limit of the general 
1-loop results of Drees and Nojiri [22] for the MSSM, but is of less interest there as it 
requires two squarks nearly degenerate with a neutralino LSP - a somewhat unusual 
spectrum. 

For completeness, it is also worth remarking that if Xi an d X2 are both fermionic as in 
models of type A, the existence of a light boson would allow the amplitude to be independent 
of m x . This follows from the same non-relativistic viewpoint as above, with the exception 
that now the heavy internal line is fermionic and so the factors of m x cancel out. In practice, 
this situation is of less interest as there are no light scalars available to provide the resonant 
diagram. However, this viewpoint does explain the corresponding results of [23l [TT] , where 
the W^-loop corrections to the scattering cross-section were obtained and are explicitly inde- 
pendent of 771dm in the limit of large dark matter mass. This scaling is entirely consistent 
with the argument above, where the role of "light boson" is played by the W. We note in 
passing that in the model of [11] one should expect a significant difference in the scattering of 
Majorana and Dirac type dark matter. In the latter case, the cross section is enhanced due 
to a significant contribution from the radiatively induced charge radius. (See also a related 
discussion in Ref. (23].) 



4 Resonant WIMP annihilation 

The pseudo-degeneracy in the WIMP sector was motivated in part by the generic need for 
enhanced annihilation mechanisms for heavy WIMPs. Although we have not discussed this 
in detail here, the present scenario with nearly degenerate charged states clearly leads to 
an interesting coannihilation scenario for which the vector-mediated cross-section is further 
enhanced by the required resummation of Coulomb effects, as discussed recently in [25J. 
Moreover, for small Am this will tend to a resonant process mediated by the formation of a 
metastable (X2X2) state. 

We will develop this possibility further in this section considering the following process 

X1 + X1 — (X2X2) ->7 + 7, (28) 
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which leads to the primary emission of monochromatic gamma rays inside e.g. the galactic 
center and can be effectively probed with existing and planned experiments. The Coulomb 
enhancements for this process in the near-threshold regime were studied previously in [26j, 
but here we will focus on the resonant contribution. 

The binding energy of the {X2X2 ) bound state is a 2 m x /4. For TeV-scale WIMPs this 
is about 13 MeV, and thus can be quite comparable with the energy excess of 2Am = 
2(m X2 — m Xl ). In what follows, we analyze the regime where the IS" state for the (X2X2) 
system is just above the energy threshold, and accessible via the kinetic energy of two 
WIMPs, 



2 

< 2Am - < m x v 2 m , (29) 

where v m is the maximal allowed WIMP velocity within the halo, on the order of 3 x 1CT 3 . 
This kinematic regime is a narrow strip on the (m x , Am) plane, 



V 10 MeV J V m x J 

The spread in WIMP kinetic energy is on the order of an MeV and thus is much wider 
than both the capture width T in and the outgoing width Y out = T2 7 , so that the Breit-Wigner 
cross section is effectively a delta-function, 

2vr 2 

a = —gJ^E - E R ), (31) 

2 

where Er = 2 Am — m * a . In Eq. (1311) we took into account that the entrance width is 
much larger than T2 7 , assuming no significant suppression of the WIMP-matter coupling, 
and therefore the rate is controlled by the smallest width. The annihilation rate into 27's of 
a bound state of two scalar particles is given by 

= 5 10 -6 MeV x ( JJbc\ . (32) 

7 4 VI TeV/ V ; 

For the parapositronium-like annihilation of two fermionic X2 S > the rate is twice larger, and 
therefore the width T2 7 has rather minimal model-dependence. 

The capture rate is also straightforwardly calculable, but in contrast is significantly 
model-dependent. For models of Type A, the scaling of the entrance width is 



r in ~ G 2 F {am x Ym x \ ^^- - 1. (33) 



ma 2 



This rate is heavily dependent on m x , but for most of the parameter space relevant for 
WIMPs is larger than (I3"2"|) . The same holds true for capture in models of Type B. It is 
interesting to note that for capture processes mediated by the Standard Model fermions, the 
P-wave capture may become significant. It is typically assumed that in the Galactic envi- 
ronment, where the velocity of colliding WIMPs is on the order of 10~ 3 , p-wave annihilation 
is suppressed relative to the s-wave by a v 2 factor, or 10~ 6 . In fact, the suppression factor 
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is far less dramatic and for models with r-exchange it is m^t^/m 2 , which is on the order of 
10~ 2 for a 100 GeV WIMP. This allows for a sizable recombination rate of a neutralino pair 
with quantum numbers L = 1; 5 = 1; J = into a bound state of two staus. 

We are now ready to estimate the total annihilation rate by averaging ( 13TI) over the dark 
matter velocity distribution. If this distribution is approximated as Maxwellian, then the 
rate is simply given by 

/ 47T \ 3/2 

(crv) = [ — g*r 2j exp(-E R /T eS ), (34) 

where T cS must be identified with the effective 'temperature' of the WIMP gas, |T eff = 
(E). For the typical velocity of dark matter inside the halo, this effective temperature is 
approximately 

T eff ~ 300 keV x (YTeV) " ( 35 ) 
Inserting some characteristic numbers for a TeV mass WIMP, we obtain 

— exp (-E R /T eS ) , (36) 
m x J 

which for Am tuned so that E R and T e fj are comparable, can be very large and indeed several 
orders of magnitude larger than the characteristic mono-energetic 7 signal [27] in generic 
MSSM scenarios for example; similar enhancements were obtained in (26] via a resummation 
of Coulomb effects near the bound-state threshold. Such mono-energetic signals are far less 
sensitive to astrophysical backgrounds and so this regime is clearly promising for indirect 
searches. Note that such a large galactic annihilation rate (1361) is not in contradiction with 
the required picobarn rate at freeze out, as this resonant effect is less important in the 
latter regime due to the larger effective temperature, which leads to an O(10 -2 ) — O(10~ 3 ) 
suppression in the thermal average. Finally, we should note that this form of enhancement 
can also affect the production of charged states in the annihilation process, although the 
rate would be more model-dependent than for monochromatic photons. 



5 Concluding Remarks 

In this paper we have embarked on a study of some of the novel direct and indirect detection 
signatures of pseudo-degenerate WIMP dark matter. Although this scenario lies within the 
characteristic WIMP sector of viable supersymmetric models and is thus far from exotic, 
the MeV-scale degeneracy in the spectrum does lead to some rather exotic signatures; most 
prominently through recombination processes with nuclei. In this concluding section, we will 
briefly summarize the results and comment on some further probes of this scenario. 

Although we have refrained from performing detailed calculations of the capture rates, 
and our treatment of the nuclear physics was somewhat cursory, the estimates obtained in 
Section 3 are sufficient to draw some broad conclusions on direct detection to focus future 
studies. For models of Type A with a weak current, we found that searches for anomalous 
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heavy isotopes in fact do not impose any significant constraints on the inelastic cross-section 
with nuclei. However, our estimate for the rate was relatively small, at best comparable 
with the existing direct detection limit on the elastic cross-section. For models of Type B, 
the cross-section is constrained by heavy isotope searches, but only in a rather small range 
arising from binding with light elements. For splittings in the range (jlj), the constraints 
are far weaker and our estimate for the capture rate suggests that simple models could 
produce rates a few orders of magnitude larger than the existing constraint on the elastic 
cross-section. Clearly it would be very interesting to explore this possibility further and to 
determine whether conventional direct detection experiments could probe capture processes 
of this type. 

A significant motivating factor in the analysis of this scenario was that the novel direct 
detection signatures may imply somewhat weaker constraints on the cross-section with mat- 
ter than would apply to generic WIMPs. This opens up the possibility of other interesting 
indirect detection possibilities in galactic astrophysics. A possible signal of this type, the 
observed excess of 511 keV photons from the galactic center, was considered recently within 
this framework |28j but seemingly required too large a cross-section to be consistent with 
terrestrial heavy isotope bounds (in this case the mass splitting required was 0(1 MeV) 
which allows for binding with C, N etc). Here we have explored several other signatures for 
indirect detection. We noted that the possibility of observing highly energetic neutrinos from 
DM annihilation in the sun is significantly diminished within these scenarios, due to the far 
higher rate for the WIMPs to be captured by iron than to annihilate. If the mass splitting 
is tuned appropriately, we also observed that the annihilation rate to two mono-energetic 
7's, e.g. in the galactic center, could be resonantly enhanced by several orders of magnitude 
through an intermediate (xt XsT ) bound state. 

It would clearly be interesting to explore other indirect signatures beyond the traditional 
7-signal from WIMP annihilation, as may arise through the possibility for recombination 
processes in astrophysical environments which may lead to new detection strategies, or al- 
ternatively more stringent constraints to be placed on these scenarios. Finally, we should 
emphasize that we have focused here on direct and indirect detection of WIMPs in the galac- 
tic halo, and have avoided the question of direct production in colliders. This is primarily 
because direct production would at least in the near term be limited to the low mass regime, 
which contrasts somewhat with part of the motivation for the pseudo-degenerate scenario, 
which was the need for heavier states. Nonetheless, this question is certainly of general in- 
terest and existing analyses of the collider physics of the WMAP strip region in the CMSSM 
should include the regime where the neutralino LSP is a pseudo-degenerate WIMP. In this 
context, the charged xf states could be sufficiently long lived to either escape or become 
bound inside the detector. 
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